Non-ideal artificial phase discontinuity in long Josephson 0-K-junctions 
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We investigate the creation of an arbitrary K-discontinuity of the Josephson phase in a long Nb- 
A10a;-Nb Josephson junction (LJJ) using a pair of tiny current injectors, and study the formation 
of fractional vortices formed at this discontinuity. The current 7inj, flowing from one injector to the 
other, creates a phase discontinuity k oc Jinj. The cahbration of injectors is discussed in detail. The 
small but finite size of injectors leads to some deviations of the properties of such a 0-k-LJJ from the 
properties of a LJJ with an ideal K-discontinuity. These experimentally observed deviations in the 
dependence of the critical current on Jinj and magnetic field can be well reproduced by numerical 
simulation assuming a finite injector size. The physical origin of these deviations is discussed. 

PACS numbers: 74.50.+r, 85.25.Cp 74.20.Rp 
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I. INTRODUCTION 

It has been shown theoretically and experimentally 
that by utilizing unconventional superconductors or fer- 
romagnetic barriers one can fabricate so-called Joseph- 
son 7r-junction3i*2i^i^(7r JJs). While for conventional 0- 
j unctions the first Josephson relation reads Is = Ic sin^, 
for TT-junctions it is Is = —Ic sin ji = Ic sin(/i -|- tt). 

In a ID long Josephson junction (LJJ) with al- 
ternating and TT regions vortices carrying half of 
a magnetic flux quantum $0 (so-called semifluxons^) 
can spontaneously appear at the boundaries between 
and TT regions Such semifluxons were observed 
experimentallyi^iiSiiiiiSii^ 

In contrast to an integer Josephson vortex (fluxon), 
which is a free moving soliton, a semifluxon is like a 
spin i particle — it is pinned at a O-tt boundary, has 
one of two possible polarities ±|^'o and, the most im- 
portant, can form the groundstate of the system. Semi- 
fluxons are very interesting non-linear objects: they can 
form a variety of groundstateSfi^4i^*iSiiLi& may flip^SiiiiiS 
rearrangai^ or get depinnediiiSfliSi by a bias current, lead 
to half-integer zero-field steps^SiSLS^ and have a charac- 
teristic eigenfrequency^ Huge arrays of fractional flux 
quantaMii^ii^^ were realized and predicted to have a tun- 
able plasmon band structure24 which can be thought of 
as a plasmonic crystal similar to photonic crystals. Semi- 
fluxons are also promising candidates for storage devices 
in classical or quantum domain. They have a very small 
effective mass^ and the potential for the observation of 
macroscopic quantum effects and building qubits. 

Recently we have demonstrated that one can also study 
semifluxons in a conventional Nb-AlO^-Nb LJJm^ This is 
done by creating a discontinuity of the Josephson phase 
(j) (the difference between and /i is explained below) 
by means of current injectors. Two injectors of width 
Aw are attached to one of the junction's electrodes at 



FIG. 1: Sketch of a Josephson junction with a pair of current 
injectors to create an arbitrary discontinuity of the Josephson 
phase. 



a distance Aa; from each other, as shown in Fig. ^ To 
create an ideal discontinuity, both the injector width Aw 
and the distance Ax between them must be much smaller 
than the Josephson penetration depth Aj. The current 
/inj, that is passed from one injector to another creates a 
phase shift k cx /jnj across the small part ~ 2 Aw + Ax of 
the top electrode. Since k cx /inj, the phase shift k can be 
tuned to K = TT by changing Imj ■ When the discontinuity 
of the phase is created, the fractional vortex centered at 
the discontinuity immediately appears to compensate it. 

LJJs with such artificially created discontinuities have 
a number of advantages. First of all, by creating a dis- 
continuity K 7^ TT one can study not only semifluxons, 
but vortices with any arbitrary flux^ The variable k 
can also be used as tuning parameter in some devices, 
6.17., to tune the plasmon band structured Second, due 
to low damping in Nb-AlO^^-Nb LJJs one can study the 
dynamics of the fractional vortices. Exponentially low 
damping at T ^ Tc due to the energy gap also helps to 
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build qubits with good dccohcrcncc figures. 

LJJs with double injector structure can also be used 
for the controllable insertion of an integer fluxon (k = 
2tt) into the LJJ^SSiSi e.g., to set/read out the state of 
a fluxon or semifluxon based (qu)bit>i2i2242S Thus, the 
structure shown in Fig. ^ can be used in a range of the 
electronic devices as well as for basic studies involving 
integer and fractional Josephson vortices. 

In experiments, one of the main steps is the calibration 
of the injectors, i.e., the experimental determination of 
the ratio «;//inj- Although a direct theoretical calculation 
of the ratio n/hni is possible in the framework of a ID 
model presented below, in practice a different approach 
is easier and more reliable: one measures the 7c(/inj) de- 
pendence and compares it with the /c(k) dependence ob- 
tained theoretically. 

The purpose of this paper is to present this calibra- 
tion method, i.e., study how the /c(k) and, therefore, 
the /c(^inj) dependence should look like for typical ex- 
perimental parameters. We especially are interested in 
the effects caused by the finite injector sizes, i.e., when 
the discontinuity is not an ideal step function. 

This paper is organized as follows. In section^ we in- 
troduce the model of a LJJ with phase discontinuities. In 
section Hill wc first consider (5-like current injectors that 
correspond to an ideal step-like discontinuity (or discon- 
tinuity point) of the Josephson phase and discuss the 
/c(k) dependence. In section llVI we compare our the- 
ory with experimental results. As there are certain mis- 
fits between measurements and theoretical predictions, in 
Scc.|3we take the finite size of the injectors into account. 
Namely, we investigate the effect of the finite injector size 
on the Ic{k) dependence and on the dependence of Ic on 
magnetic field H and compare them with experimental 
results. Section IvH concludes this work. 
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FIG. 2: Functions 9{x), 9x{x), 9xx{x) and the approximation 
79(2;) of 9xx{x) by two rectangular pulses of width Aw and 
distance Aa;. 



It is clear from Eq. ^ that the solution <f>{x) is also 
TT-discontinuous at the same points as 6{x). Therefore, 
we often call the points where 0- and 7r-parts join phase 
discontinuity points. 

The bias current density 7(2;) and 0xx{x) are two ad- 
ditive terms in Eq. This means that if one does 
not have initially a 9xx{x) term (like in a conventional 
LJJ), its effect may be substituted by an additional bias 
current je{x) = Oxx{x). If 0{x) = 7rH(x) is a step func- 
tion [Fig. |2a)], where H(x) is the Heaviside step func- 
tion, then Ox{x) = ■n5{x) is a (5- function of strength tt 
[Fig.I3;b)], and 9xx{x) = -n5{x)/x [Fig. I^^c)]. Thus the 
TT-discontinuity of the phase corresponds to (or can be 
substituted by) an additional current with the density 
~TTd{x)/x, i.e., by a current dipole of strength Tj-^SiSSi^ 
Of course, in experiment one cannot pass infinitely large 
currents via infinitesimal wires, so the real current pro- 
file will be an approximation of Oxx{x), e.g., like the 
one shown in Fig.|2Id). Nevertheless, following Ref . WK. 
we first consider an ideal current dipole ((5-like current 
dipole) corresponding to an ideal discontinuity. 



II. THE MODEL 

The dynamics of the Josephson phase in a LJJ consist- 
ing of 0- and 7r-parts can be described by the ID time- 
dependent perturbed sine-Gordon cquationSiSSiSi 

fpxx - 4>tt - sin0 = a(j)t - 7(x) - 9xx{x), (1) 

where 4'{x, t) is the Josephson phase and subscripts x and 
t denote the derivatives with respect to coordinate x and 
time t. In Eq. the spatial coordinate x is normalized 
to the Josephson penetration depth Xj and the time t 
is normalized to the inverse Josephson plasma frequency 
u!~^; a = l/^/JJ^ is the dimensionless damping parame- 
ter; 7 ~ j / ic is the external bias current density normal- 
ized to the critical current density of the junction. The 
function 9{x) is a step function which is discontinuous 
at the points where the 0- and 7r-parts join and constant 
elsewhere. For example in O-tt-LJJs, 0{x) can be equal 
to zero along all 0-parts and to tt along all 7r-parts. For 
simplicity we assume that there is only one discontinuity 
point at the center of the junction, a; = 0. 



It is clear, that when we speak in the language of cur- 
rent dipoles wc can change the amplitude of the cur- 
rent, i.e., use ^k,S{x)/x, which corresponds to a n- 
discontinuity of the Josephson phase (instead of a tt- 
discontinuity). Below wc consider a LJJ with such an 
arbitrary ^-discontinuity of the Josephson phase, i.e., 
0{x) = Kn{x) in Eq. ||T}. 

Note that other authora^iiSiSS often use directly the equa- 
tion written for the continuous phase fi{x,t), which for 
the case of a ^-discontinuity reads 



t^xx - Mtt sin in + 9{x)) = a/it - l{x). (2) 



Equation and |(2J) are, actually, equivalent and 
one can be obtained from the other by substitution 
4>{x,t) = fJ-{x,t) + 9{x)& To simulate the system we 
used the "StkJJ" software^Si which implements finite 
difference algorithms for solving Eq. Q. 
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III. IDEAL DISCONTINUITY 

Using Eq. we first would like to know how the nor- 
malized maximum supercurrent 

/c(0,0) 

depends on k. Here v can be obtained by varying the 
normalized supercurrent 

1 /-'/^ 
i = J sincj) dx 

« J-l/2 

with respect to ^, which should still be a solution of 
Eq. QJ. The length / = L/Xj is the normalized LJJ 
length. So far the idn) dependence was calculated only 
for an annular LJJi2L2S It was predicted and measured2& 
that this dependence is a Fraunhofer pattern with the 
first minimum at k = 27r. This result does not depend on 
the junction length L. 

In a junction of linear geometry the situation is differ- 
ent and Jc(t) depends on the normalized length I. Let us 
first consider extreme cases. 

1. / ^ 1. In this case one can follow the derivation 
of Ref. and arrive to the same Fraunhofer pattern (cf. 
Eq. (10) of Ref. m 



sin(K/2) 



k/2 



(3) 



On the other hand it is clear, that a ^-discontinuity and a 
(K-|-27r)-discontinuity in linear LJJs are physically equiv- 
alent, since the phase is defined modulo 27r and there are 
no topological limits as in annular LJJs. Therefore, the 
Fraunhofer pattern (PJ should be replicated periodically 
with a period of 2tt along the K-axis, as shown in Fig. |2| 
(dotted lines). Hence, for a given value of n one may 
have more than one ic, corresponding to different solu- 
tions. If one sweeps k at fixed i, or sweeps i at fixed 
K the system may switch from one solution to another 
by emitting a fluxon which leaves the junction and does 
not affect the system later. This happens when the bias 
point crosses one of the dotted lines in Fig. 13 We note 
that in an annular LJJ a fiuxon cannot escape. Therefore 
transitions as described above (for a linear LJJ) cannot 
occur, and only the principal Fraunhofer curve Eq. Q 
has to be considered. We do not discuss the stability of 
each solution here, but it is important to note that the 
topmost branch, indicated by the solid line in Fig.|2| sep- 
arates all Meissner states from the finite voltage state and 
can be measured in experiment. Thus, one gets cusp-like 
minima at k = t: + 2Tm (with integer n) . 

2. Z <C 1. In case of a short junction one can neglect 
the spatial variation of the Josephson phase along the 
junction and assume that the phase in the left half of 
the junction is (po and in the right half of the junction is 
(/)o -l- K. By varying (f>o we find that the critical current 




n 2k 3jt 

phase discontinuity k 



4jt 



FIG. 3: Dependencies of tlie normalized critical current ic 
on the value of the Josephson phase discontinuity k. Dotted 
lines show Fraunhofer patterns shifted by 0, 27r, 47r, the 
corresponding ic(K) for L = 00 is shown by solid lines. Dashed 
lines show the dependence @ for L Aj. 



depends on k as 



(4) 



which is shown in Fig. O by the dashed line. 

In Fig. one can see that both dependences have 
a common feature: they have maxima at k = 27rn 
and cusp-like minima (possibly with hysteresis) at k = 
TT + 2'nn. The value of v at the minima depends on the 
junction length and varies from for a short junction to 
2/7r for an infinitely long one. 

3. Intermediate /. In this case the ic{i^) curve lies in 
between the two curves for I <C 1 and i ^ 1, but still 
has cusp- like minima at k = 7r(2n -1-1). The numerically 
calculated dependence of the minimum critical current 
ic(7r) on junction length will be shown later in Sec. 
Since cusp- like minima in ic{n) always correspond to a 
phase discontinuity k = (2?i+l)7r, the calibration process 
of the injectors is simple. One should measure /c(/inj)- 
Since k oc /inj the value of /inj at the first cusp-like min- 
imum of the /c(Anj) dependence corresponds to k mod 

2-K = -K . 



IV. SAMPLES 

So far an ideal K-discontinuity or, to put it in terms of 
an additional bias current, (5-like current injectors were 
discussed and the ic{n) dependence was investigated. As 
a result we saw that the ic{n) curve is 2-k periodic and 
has a minimum of the critical current for phase disconti- 
nuities K = {2n + l)7r independent of the junction length. 
In order to experimentally investigate artificial 0-k-LJJs 
we used samples that were fabricated at the University 
of Karlsruhe as well as by HypreaS and that are based 
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FIG. 4: Experimentally measured /c(J'inj) dependence at H = 
(sample #2). 

on conventional Nb-AlO^^-Nb technology. All data pre- 
sented here have been obtained at T = 4.2 K from dif- 
ferent samples with parameters summarized in Tab. 
For all samples the junction layout is identical to the one 
presented in Fig.^ A photograph can be found in Fig. 3 
of Ref. [2^ (sample #4 in Tab.Hl). However, the values of 
Aj, the lengths L and the injector sizes are different. For 
all samples the I-V characteristics (IVC) and the depen- 
dences of Ic{H) without injector current were measured 
to ensure good sample quality and the absence of trapped 
magnetic flux (not shown). 

Figure 0] shows a typical measurement of Ic (/inj ) at 
H = 0, obtained from sample #2. One can clearly see 
the almost periodic modulation of the critical current as 
a function of /jnj which is very similar to the one pre- 
dicted in Fig. 13 Yet, there is an unexpected decrease of 
Ic at the side maxima /inj ~ ±19 mA. We have also mea- 
sured the Ic [H) dependence at /mj corresponding to the 
first minimum in the /c(/inj) dependence, i.e., the Ic{H) 
dependence in the 0-7r-state. Such a dependence is ex- 
pected to have a minimum of Ic at H = 0ii2i2^ In Fig. [S] 
two typical Ic{H) dependences experimentally obtained 
from different samples are presented. Although we in- 
deed observe a minimum of Ic at H = 0(^*^ the curves 
in Fig.Elare rather asymmetric with respect to H = 0*2^ 
the amplitudes of the first maxima are different and the 
higher order maxima have different, sometimes not very 
regular, periods. 

In the following we show that these deviations can 
be reproduced in numerical simulations using the model 
with the finite size of the injectors. The physical origin 
of these deviations is discussed. 



V. NON-IDEAL DISCONTINUITY 

To model the finite-sized injectors we used the cur- 
rent injection profile 7e(a;) shown in Fig.|2Id). Note that 



this shape is a rather simple, but decent approximation, 
corresponding to our system with superconducting injec- 
tors and a ground plane. Also, we consider here a ID 
model, whereas in the real geometry (see Fig. ^ or Fig. 
3 of Ref. |23) the injector current may be distributed in a 
rather peculiar way across the top electrode. 
Equation ^ can be rewritten as 

(j)xx- (i)tt-iivcL4'^ acj)t~"f -^e{x), (5) 

where ^e{x) approximates Oxx{x). By introducing a finite 
injector size the K-discontinuity is now stretched over a 
finite distance = 2Azi; -I- Ax. In the injector area 
the dynamics of the Josephson phase and therefore the 
junction properties are dominated by the injector current 
distribution. 

First, we are interested in the static phase profile in the 
injector area. If the injectors provide the current density 
distribution je{x), the static sinc-Gordon equation is 

(l)xx = sine/) - J - j0{x). (6) 

For the following considerations it is more convenient to 
separate the amplitude 7inj of the injector current density 
from the actual injection profile 

r I , -f <x< 

le{x) = 1 , < X < +^ , (7) 

[-1 , +^ <x< +§ 

so that 76/ (x) = 7inj7e(2^)- In the injector region, 7e(x) 
provides the dominating contribution (|7inj| 3> 1) to the 
rhs. of Eq. ©. If in addition Ax is assumed to be <C 1, 
Eq. lO can be simplified to 

4>xx = -7inj7e(a;)- (8) 

Double integration of Eq. © gives the phase profile 0(.t) 
in the region of injectors (— A^/2 < x < -I-A//2) depicted 
in Fig.El The difference of the phase cj){Al/2)-(j){-Al/2) 
is proportional to 7inj and should be equal to k. Thus 
we obtain a relation between k and the amplitude of the 
injector current density 7inj 

K = -7i„jAw(Aw + Ax). (9) 

Now, let us consider the effect of a finite injector size 
on the ici^) dependence. Fig.Qshows two sets of numer- 
ically calculated Zc(k) curves. In Fig. |7{a) junctions of 
the same length L = 2Aj but different injector sizes are 
compared, ranging from Ax = 0, Aw = (ideal discon- 
tinuity) to Aw = 0.32, Ax = 0.33. The most prominent 
difference compared to an ideal discontinuity is a decrease 
in critical current at the side maxima ic{2Tm) {n ^ 0). 
The difference of critical currents ic(0) — ic(27i7r) becomes 
smaller as the injector size is reduced and approaches zero 
as injectors become ideal (Ax 0, Aw ^0). In addition 
the minima of the curve shift slightly towards values of 
|k| > 7r(2n-|- 1) and ic,min is slightly reduced when injec- 
tors become very large and comparable to the Josephson 
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"fabricated at the University of Karlsruhe, Germany, 
''fabricated by Hypres, Elmsford, New York, USA. 

TABLE I: Sample properties. The quantities I, Aw and Ax are the normalized (with respect to Aj) LJJ length, injector width 
and distance between the injectors. N is the number of junctions on the chip, that were measured. The junction width is 5 fim 
for all samples. 
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FIG. 5: Measured Ic{H) dependences in the O-vr-state for two LJJ (chip #1) of length L = 1.25 Aj (a) and L — 3.75 Aj (b). 
Arrows mark the crossing of the _ff-axis in negative fields [cf. Fig.[Tl(b)]. 
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FIG. 6: Spatial dependence of 7e(a;) (solid line) and 4i{x) 
(dotted line) in the injector area. Aa; = Aw — 1, 7inj = —n/2 
for simplicity. The injector size is intentionally chosen to be 
huge (~ Aj) to visualize the behavior of the Josephson phase 
in all regions. 

penetration depth. The last effect is rather small and 
can normally be neglected for typical sizes of injectors. 
Fig. Elb) shows the id^) dependence of junctions with 
the same injectors, but for different LJJ lengths. Again, 
there is a decrease of ic(2mT) for n ^ 0; yet, with increas- 



ing junction length L the ic{K) dependence approaches 
the curve of an ideal discontinuity (sec Fig.|2|). 

Fig. |S1 shows the numerically calculated dependence of 
the minimum critical current ic,min at k « tt on LJJ 
length L, for ideal and non-ideal discontinuities, as well 
as experimentally obtained results. The numerical calcu- 
lations for the case of finite injector size were performed 
with Aw = 0.2, Aa; = 0.21, comparable to the experi- 
mental dimensions of sample series ^1. In the simula- 
tions the finite injector size has a negligible influence on 
*c,min- However, for L < lOAj the experimentally deter- 
mined values of ic,min are systematically smaller than 
theory predicts. Only for very long JJs, L > 10 Aj, 
measurements and theory approach the same asymptotic 
value of 2/7r (semifluxon dcpinning currenti^iiSiM) . It is 
interesting to note, that the measured values (except for 
sample ^^2, that has smaller injectors) show a similar de- 
pendence as in theory, yet the experimental data seem to 
be shifted along the L-axis, towards longer JJs. So far 
this behavior is not fully understood. However the bet- 
ter agreement between theory and experimental data for 
sample #2 gives already an indication, that the actual 
injector size does play an important role and that the ID 
model for the injectors could be too simple. 
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FIG. 7: Normalized critical current ic vs. the Josephson 
phase discontinuity k, numerically calculated for junctions 
with fixed L — 2Aj and different injectors (a) and junctions 
with same injectors (Aw = 0.08, Ax = 0.09) but different 
L (b). In (a) the sizes of injectors are Ax = 0,Aw — 
(gray circles), Aw = 0.04, Ax = 0.05 (solid line). Aw = 0.16, 
Aa; = 0.17 (open circles) and Aw = 0.32, Ax = 0.33 (black 
circles) . 



Looking at the maxima in Fig. [T] again, we conclude 
that the deviations from the ideal characteristics depend 
on the relative size of the injector region with respect 
to the junction length. There is a simple picture to 
understand this qualitatively. For example, at k ~ 2tt 
the Josephson phase increases by 27r, from (j>{~Al/2) to 
(j){+Al/2) in the region between the injectors, so that in 
this region the net Josephson current across the junction 
can be neglected (if the increase in (f> is linear) , no matter 
what bias current is applied, because 7 ^ 7inj. There- 
fore the junction length L is effectively reduced by the 
size of the injector region and to a first order approxima- 
tion the critical current is ic(±27r) w 1 — (Aw + Ax) /I. 
Note, that the effective injector region, that has a net 
Josephson current equal to zero, is only Ax + Aw rather 
than Ax + 2Aw, due to the parabolic shape of 4> inside 
the injector electrodes (cf. Fig. IBJ. In Fig. |21 we plot 
Aic = 1 — jc(27r) vs. {Aw + Ax) /I obtained experimen- 
tally and numerically for junctions with different injector 
sizes and of different LJJ lengths I. One can notice, that 
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ideal discontinuity 
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FIG. 8; Minimum normalized critical current ic vs. junction 
length numerically calculated for an ideal discontinuity (solid 
line) and non-ideal injectors with Ax — 0.21, Aw — 0.2 (solid 
circles). Open symbols correspond to measurements of the 
samples summarized in Tab.|T| 
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FIG. 9: Numerically calculated (filled circles) and measured 
(open symbols) reduction of normalized critical current Aic = 
1 — ic(27r) for different junction lengths and injector sizes. 



the experimental data are in a good agreement with nu- 
merical calculations, yet the measured drop in critical 
current is always slightly larger than the theory predicts. 
By confining even larger phase jumps (k = 47r, Gtt, . . .) 
into the area of injectors, the critical current of the junc- 
tion slowly decreases further (sec Fig. [TJ. Yet, the rela- 
tive change ic{2Trn) — ic(27r(n-|- 1)), n > is much smaller 
than Aic = 1 — *c(±27r) and depends on the details of the 
phase bending inside the injector area. This, however is 
a higher order effect which is not addressed in this pa- 
per, because in experiments one always wants to have 
Ax, Aw < 1. 

Now let us concentrate on the 0-7r-state [/jnj corre- 
sponding to a minimum in /c(-^inj)] and consider the 
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FIG. 10: Numerically calculated ic{h) dependence of a LJJ 
[L — 2Aj) in the 0-7r-state. Injector size is Aw = 0.16 and 
Aa; = 0.13. 



icih) dependence, where h = 2H/Hci and Hd = 
$o/(27r/ioALAj) is the first critical field for a LJJ with 
thick electrodes. In Fig. ^| we show the numerically cal- 
culated ic{h) dependence of an artificial O-tt-LJJ with 
finite injectors. Parameters arc chosen to be compara- 
ble with experimental dimensions^ (cf. sample #4). As 
has already been mentioned in Ref. |22, simulated results 
reproduce all features of the experimentally measured 
Ic{H) dependence rather well [see Fig. El^b) or Fig. 5 
in Ref. |22|. Although the curve shown in Fig. IIUI has 
the characteristic minimum of the critical current in zero 
magnetic field, it is rather asymmetric. The maximum 
critical currents of the side lobes arc different and the pe- 
riodicities of the outer oscillations for negative and pos- 
itive field differ. Yet, with decreasing injector size this 
asymmetry disappears as can be seen in Fig. 1111 

To understand the reason for this asymmetries, we, for 
the moment, consider a relatively short JJ (L < Aj) and 
find the expression for ic{h) for injectors having a nor- 
malized size Aw and distance Ax between them. The 
Josephson phase is assumed to be linear outside the in- 
jector area, i.e., 



4>a + h{x ■ 
bo + K. + hix 
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Inside the injector region the phase changes from to 
(^0 + over the distance A/ — 2Au' 4- Ax, but its partic- 
ular profile will be discussed later. 

The total supercurrent carried by the LJJ outside the 
injector area is calculated by integrating sin^o„t(a;) over 
the junction length excluding the injector region. By 
varying i^g to maximize we find that the maximum 
is reached at 
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FIG. 11: Numerically calculated dependence of a LJJ 

(L — 2Aj) in the 0-7r-state. Injectors have a width/distance of 
Aw = O.Ol/Ao; = 0.011 (open circles). Aw = 0.04/As: = 0.05 
(black filled circles) and Ato = 0.08/Aa: = 0.09 (gray filled 
circles) . 



and is given by 



= _ cos - 



h{l-Al) . h{l-Al) 
— sin . 



(12) 



The dependence (|12f) is an odd function of h, i.e., it is 
anti-symmetric with respect to the origin h = 0. We 
note that in experiment one always measures the absolute 
value |ic(/i)|, which is a symmetric (even) function with 
respect to the origin h = 0, and has a minimum at /i = 0. 

Now let us calculate the supercurrent carried by the 
region of injectors. The phase 4>{x) changes from (j)Q to 
(/)o + K as shown in Fig. El i.e., it has two parabolic parts 
corresponding to two injectors and a linear part corre- 
sponding to the spacing between injectors. To calculate 
the supercurrent is.inj carried by the injector region, we 
again integrate 



1 



-Al/2 



sindlix) dx. 



Al/2 



(13) 



Since some parts of (f>{x) are parabolic, we will end up 
with Fresnel integrals which will only add additional com- 
plications. Wc thus change the phase profile a little bit 
to simplify integration, taking into account that the in- 
jector current profile shown in Fig.|2]is only a simple ap- 
proximation. Considering that in our original model (t>{x) 
changes parabolically, i.e., slowly, and that for supercon- 
ducting injectors most of the current flows close to the 
inner edge of the injector, wc adopt a very simple injector 
phase profile 
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/ X l\ Ax + Aw Aa; + 



I.e., the phase changes hnearly from (pQ to </)o + k only 
in the space between the injector centers. By effectively 
shortening the injector area, the still unattended outer 
halves of the injectors are simply added to the part of the 
LJJ without injectors, i.e., we substitute A/ = Aw + Ax 
for A/ in Eqs. and (USJ. Thus, the supercurrent 

carried by the injector region calculated for optimum 
00 (dU is 

2 K 

is.inj = (Ax + A-w)— sin-. (15) 

Note, that is,inj does not depend on magnetic field, be- 
cause we did not take the penetration of magnetic field 
into the injector region into account. 

The resulting normalized maximum supercurrent vs. 
field dependence which one measures in experiment is 
given by 

is[h) = \ic{h) + is,i„j|, (16) 

i.e., the anti-symmetric ic{h) dependence from Eq. (|12|l is 
first shifted along the i-axis by is,inj and then the absolute 
value is calculated, [see Fig.^](a)]. As a result one gets 
asymmetric maxima for h > Q and /i < 0. Moreover, the 
side oscillations of ic{h) + is.inj due to the shift may not 
intersect the /i-axis anymore for ft, < 0, but do intersect 
it for ft > 0, resulting in very different patterns of the 
side maxima for the i,,(ft) dependence, [see Fig.^Ja)]. 

Before we proceed with a comparison of experiment 
with our analytical results, we would like to note an- 
other effect which may become noticeable in not very 
long junctions with rather large injectors. Namely, the 
phase jump n which, up to now, was supposed to be only 
due to the injector current k = Kjnj, can also have a 
contribution from the magnetic field penetrating the in- 
jector area. We write k = Kinj + ft (Ax -I- Aio) and use 
the so-defined k in Eq. (|15|l and (|l(j|) . Then, the contri- 
bution of the injector area to the supercurrent, i.e., the 
above mentioned shift, becomes field dependent. This 
means that side oscillations of Zc(ft) + is,inj which did not 
cross the {ig = 0)-axis for lower fields, may do so for 
larger fields, as is shown in Fig. I12r b'). The larger the 
ratio (Ax-I- Aw)/^ the more noticeable this effect will be. 
Fig.Elshows experimentally obtained /c(^)-depcndences 
of two LJJ of different length in the 0-7r-state. As one 
can see clearly, these measurements are in good agree- 
ment with the above calculations. In particular, one can 
notice that the shorter the JJ is, the more pronounced is 
the asymmetry in Zc(ft). 

It is interesting to note, that with the simple short 
junction model presented above, almost all features of 



I 

the measured Ic{H) dependence can be explained and 
effectively pinpointed to the supercurrent inside the in- 
jector area. 



VI. CONCLUSIONS 

A pair of tiny current injectors can be used to cre- 
ate an arbitrary phase jump of the Josephson phase and 
to study arbitrary fractional vortices pinned at it. An 
important step, which one should make following this 
approach, is to calibrate the injectors, i.e., to measure 
the Ic{Iini) dependence. Here we calculated how /c(/inj) 
should look like in the 0-k-LJJ of length L with an ideal 
discontinuity. We discuss the effect of finite injector size 
on the /c(-finj) and Ic{H) dependence and find rather 
good agreement between our model and experimental 
data. The consequences of finite sized injectors compared 
to an ideal discontinuity can be summarized as follows: 

(I) The ic(K) dependence is not perfectly 2-k periodic 
but exhibits a decrease in ic(±2n7r) < 1, n 7^ 0. The 
drop of ic at K = ±27r can roughly be estimated as 
1 — ic (±27r) w (Aw -I- Ax) /I and is associated with the 
exclusion of injector area that carries no net supercur- 
rent. 

(II) The ic(ft) characteristic for k = tt is asymmetric 
with respect to magnetic field. This asymmetry is caused 
by the Josephson currents flowing in the area between the 
injectors. These Josephson currents shift the ic(ft) curve 
by is, inj [cf- Eq. ((T3|l ] and become weakly field dependent 
for A^ < I. As in experiment only is(ft) = |ic(ft) + is, inj | 
can be measured, yet ic{h) is an odd function of ft, parts 
of the ic(ft) curve get reflected at the ft-axis, resulting in 
a field dependent and asymmetric pattern of the higher 
order maxima. The difference in normalized critical cur- 
rent between the main maxima is « ±^(Ax+Aw)// with 
respect to ideal injectors. 

In general, it is desirable to make the injector region as 
small as possible, to avoid the discussed finite size effects. 
However, for cases where this is not possible, our results 
allow one to interpret the experimental data properly and 
optimize the design of such LJJs. 
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FIG. 12; (a): Calculated dependence of js(/i) (black line) and |i3(/i)| (gray line), (b): Calculated dependence of |is(/i)| in case 
magnetic field is penetrating the injector area. Parameters are I = 2, ft = n and Aa; = Aw = 0.2. 
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